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COMPUTING HERMITIAN DETERMINANTAL 
REPRESENTATIONS OF HYPERBOLIC CURVES 


DANIEL PLAUMANN, RAINER SINN, DAVID E SPEYER, AND CYNTHIA VINZANT 

Abstract. Every real hyperbolic form in three variables can be realized as the 
determinant of a linear net of Hermitian matrices containing a positive definite ma¬ 
trix. Such representations are an algebraic certificate for the hyperbolicity of the 
polynomial and their existence has been proved in several different ways. However, 
the resulting algorithms for computing determinantal representations are compu¬ 
tationally intensive. In this note, we present an algorithm that reduces a large 
part of the problem to linear algebra and discuss its numerical implementation. 


Introduction 

Let / be a real homogeneous polynomial of degree d in variables x, y, z. A Her¬ 
mitian determinantal representation of / is an expression 

(1) / = det(xMi + yM 2 + zM 3 ), 

where Mi, M 2 , M 3 are Hermitian d x d matrices. The representation is definite if 
there is a point e G M 3 for which the matrix ei Mi + e 2 M 2 + e 3 M 3 is positive definite. 

This imposes an immediate condition on the projective curve Vc(/)- Because the 
eigenvalues of a Hermitian matrix are real, every real line passing through e meets 
this hypersurface in only real points. A polynomial with this property is called 
hyperbolic (with respect to e). Hyperbolicity is reflected in the topology of the 
real points Vr(/). When the curve Vc(/) is smooth, / is hyperbolic if and only if 
Vr(/) consists of |_f J nested ovals, and a pseudo-line if d is odd. 




FIGURE 1. Cubic and quartic hyperbolic curves in P 2 (M). 

The Helton-Vinnikov theorem [3] (previously known as the Lax conjecture) says 
that every hyperbolic polynomial in three variables possesses a definite determi¬ 
nantal representation ([TJ) with real symmetric matrices. Thus given a hyperbolic 
plane curve, one can investigate the problem of computing a definite determinantal 
representation. 

Computing symmetric determinantal representations of hyperbolic plane curves 
both symbolically and numerically was investigated by Sturmfels and two of the 
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current authors in [5] and in the case of quartic curves in [4]. Recently, it was 
discovered [B] that looking for Hermitian matrices M \, M 2 , M- il rather than real 
symmetric matrices, greatly simplifies this computational problem and the proof of 
Helton and Vinnikov’s theorem. 

The goal of this paper is to present an algorithm for computing determinantal 
representations ([TJ) , examine this algorithm numerically, and to compare it with 
existing methods. This construction is based heavily on [B], which generalizes a 
classical construction due to Dixon |TJ. 

Acknowledgements. We would like to thank Arno Fehm for helpful discussions 
concerning Hilbert’s Irreducibility Theorem. Cynthia Vinzant was partially sup¬ 
ported by an NSF grant DMS-1204447. 


1. The Algorithm 


The input to our algorithm is a polynomial / G R[x,?/, z\ of degree d with smooth 
complex variety Vc(/) and a point e = (ei,e 2 ,e 3 ) G M 3 with respect to which / is 
hyperbolic. We will intersect the curve Vc(/) with the degree (d — 1) curve given 
by the directional derivative 
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We assume that the intersection Vc(/) H Vc (g) in P 2 is transverse. That is, the two 
curves Vc(/) and Vc(fi') intersect in d(d— 1 ) distinct points. In fact, this implies none 
of these intersection points are real [ 6 ] Lemma 2.4], If the intersection Vc (/) fl Vc(g) 
is not transverse, we may replace g with the directional derivative of / in direction 
e 7 , taking e! to be a generic point sufficiently close to e. 

The output of the algorithm will be three Hermitian dx d matrices M \, M 2l M :i 
such that / = c • det(xMi + yM 2 + zM 3 ) where c G M>o and e\M\ + e 2 M 2 + e 3 M 3 is 
positive definite. Furthermore, g will be one of the diagonal minors of the resulting 
matrix M = xM\ + yM 2 + 2 M 3 , namely the minor of M obtained by removing the 
first row and first column from M. 

The construction below is based on the idea that if the Hermitian matrix M is a 
determinantal representation of / = det(M), then its adjugate matrix M ad -> satisfies 


M adj ■ M = det(M) ■ I = f-I. 


Let a denote the top row of M ad L Then, taking the top row of this matrix equation, 
we obtain the relation a{xM\ + yM 2 + zM 3 ) = (/, 0,... ,0). Similar arguments 
give (xM\ + yM 2 + zM 3 )a T = (/, 0,... ,0) T . We introduce a suitable vector a = 
(an, 012 ,..., aid) and solve these linear equations in the entries of the Mj. This finds 
(Mi, M 2 , M 3 ) without ever explicitly computing M ad L 
The algorithm proceeds as follows: 


(Al) Compute the d(d — 1) points Vc(/) fl Vc(<?). 

(A2) Split the points into two disjoint, conjugate sets Vc(/,<?) = S U S. 

(A3) Let an equal g. 

(A4) Extend a n to a basis a = (an,..., a lrf ) of the vector space of polynomials in 
C[x,y, z]d~ 1 that vanish on the points S. (In particular, this space is always 
of dimension d; see Remark [7] below.) 
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(A5) In the 3d 2 variables (Mi)jj, ( M 3 )ij , solve the 2d(f i 1 2 2 ) = (d+2)(d+l)d 

affine linear equations coming from the polynomial vector equations 

a (xM l + yM 2 + zM 3 ) = (/, 0 ... 0) 

(xM 1 + yM 2 + zM 3 )a T = (/, 0...0) T 
(A6) Output the unique solution Mi, M 2 , M 3 . 

We need to argue that such a solution M = xM 3 + yM 2 + zM 3 exists, is unique, and 
has the desired properties, which we do below. Numerical implementation of this 
algorithm and surrounding computational issues will be discussed in Section [2] 

Theorem 1. Let f e M[a;, y, z] be hyperbolic with respect to a point e € R 3 with 
/(e) > 0. Suppose that Vc(/) is smooth and that all the intersection points ofVc{f) 
andVc(g) are transverse, where g = e i§f+ e 2 §^+ e 3 §f • Then the system of equations 
in (A5) has a unique solution M 1? M 2l M 3 , which are Hermitian matrices satisfying 

f = c ■ det(xMi + yM 2 + zM 3 ) and e^Mi + e 2 M 2 + e 3 M 3 y 0, 
where c G M >0 . 

Proof. ( Existence .) First, let us show that the affine linear equations (A5) have some 
solution Mi, M 2l M 3 . By Construction 4.5 and Theorem 4.6 of [6], there exists a Her¬ 
mitian linear matrix M' = xM[ + yM' 2 + zM 3 such that for some c / 0, the determi¬ 
nant det(M') equals c d ~ 1 f, the matrix M'(e) = eiM[-\-e 2 M! 2 -\-e 3 M 3 is either positive 
or negative definite, and the first row of the adjugate matrix A = (l/c d_2 )(M / ) ad j is 
precisely a = (an,ai 2 , • • •, aid). Since the matrices M' and (M / ) ad -’ are Hermitian, 
it follows that the first column of (l/c d_2 )(M / ) adj is a 1 = (all, an, ... ,aid) T . 

In fact, the constant c must be positive. We can see this from examining our 
matrices at the point e. Since M'(e) is definite, both (M') ad j(e) and A(e) must be 
definite as well. Furthermore, because the (1,1) entry of A(e), namely an(e) = g(e), 
is positive we see that the matrix A(e) is positive dehnite. Then the equation 

det (A) = (l/c d_2 ) d det((M / ) adj ) = (l/c d ~ 2 ) d ■ (c d ~ 1 f) d ~ 1 =c-f d - 1 

evaluated in the point e shows that c is positive. To fold a solution to the equations 
(A5), let M = (l/c)M'. Then 

det(M) = (l/c) d det(M') = (1 /c)f. 

Furthermore M adj equals (l/c) d ^ 1 (M') ad j, which is (1/c) • A. We know that both 
M adj • M and M • M ad j equal det(M)/. Dividing these identities by (1/c) we see 
that A ■ M = f ■ I and M ■ A = f ■ I. From taking the first row of the first equation 
and the first column of the second equation, we see that M satisfies the equations 

a M = (/, 0 ... 0) and M a T = (/, 0 ... 0) T . 

Since c • / = det(M) and M(e) is positive dehnite, in order to finish the proof, it 
suffices to show that this is the unique solution to these equations. 

( Uniqueness.) First we argue that if the system of equations (A5) has two distinct 
solutions, then it has two distinct solutions M = xM x + yM 2 + zM 3 and N = 
xN\ + yN 2 + zN 3 for which both det(M) and det(AI) are not identically zero. 

By the proof above, there is a solution M to the equations (A5) with det(M) ^ 0. 
If there is another solution N, then for any A € M, the matrix AM + (1 — X)N is 
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also a solution. The set of A G M for which the determinant det(AM + (1 — A )N) is 
identically zero is closed and does not contain A = 1. Thus for almost all choices of 
A 6 I, det(A M + (1 — A )N) is nonzero. 

Now suppose M = xMi + yM 2 + zM 3 and N = xNi + yN 2 + zN 3 are matrices 
satisfying the equations (A5) for which det(M) and det(iV) are not identically zero. 
We see that at a general point point ( x,y,z ) in V<c(/), the vector a is non-zero, 
and the matrix M does not have full rank. Since det(M) has degree d and / is 
irreducible, we can conclude that det(M) = af for some constant q^O. Similarly 
det(iV) = Pf for some p / 0. 

Again we use the identity (M) ad J • M — M ■ (M) ad J = det (M)I = afl. Therefore 
(M ad i)] ■ M = a(f, 0... 0) = a(a ■ M), where (M ad i)i is the first row of M ad h For 
generic (x, y, z) in C 3 , the matrix M is invertible, which implies that the first row of 
M ad ’ is aa. Similarly, the first column of M ad is ad T , and the hrst row and column 
of iV adj are Pa and Pa T , respectively. 

The matrices (l/o;)M ad j and (l/P)N ad i therefore have the same first row and 
column and both have rank-one along the curve Vc(/)■ For generic (. x,y,z ) G 
V<c(/), the entries of these rows and columns are non-zero. Therefore the difference 
(l/a)M ad j — ( y l/P)N ad i vanishes on Vc(/) and all of the entries of this difference must 
be divisible by /. However, the entries of these matrices have degree d — 1 whereas / 
has degree d, so we see that (l/a)M ad j must equal ( l/P)N ad K In particular, M ad j is 
a constant multiple of A^ ad h It follows that M is a constant multiple of N. Because 
our affine linear equations (A5) are not homogenous, we see that in fact M — N 
and our solution is unique. □ 

Remark 2. In fact, in our algorithm we can replace g in (|2]) by any polynomial 
g G M[x, y, z]d-\ where g(e) > 0 and g interlaces / with respect to e. By this, 
we mean that for every point p G M 3 the roots of the univariate polynomial g(te + p) 
interlace those of f(te +p). 

Remark 3. When the curve V(/) is singular, some of the above construction goes 
through, but requires more care. Because / has a definite determinantal represen¬ 
tation / = det(M), e.g. [6J Cor. 4.9], we know that there exists an interlacer 
g = (M ad j) n and a basis of polynomials a = ((M ad j) n ,..., (M ad j) ld ) for which 
the equations (A5) have a solution. Finding such polynomials remains a challenge. 
When / is singular, the intersection V(f) flV(g) contains points with multiplicities. 
It is unclear exactly how to split these points with multiplicity and find the correct 
linear space of polynomials an,..., that vanish on them. One example of this 
procedure is given below. 

Example 4. To illustrate our algorithm, we apply it to the quartic 

(3) f(x,y,z) = x 4 — Ax 2 y 2 T y 4 — kx 2 z 2 — 2y 2 z 2 T z 4 , 

which is hyperbolic with respect to the point e = [1 : 0 : 0], and appears as Ex¬ 
ample 4.12 in [6]. This curve has two nodes, [0:1:1] and [0 : — 1 : 1], but done 
carefully, our algorithm still works. Figure [2] shows the real curves Vm(f) and Vr(< 7 ), 
where g = §y, in different affine planes. 

First we define an to be the directional derivative \D e f = x 3 — 2 xy 2 — 2 xz 2 . The 
intersection of / and an consists of the eight points [2 : ± 1 /3 : ±i], [2 : ±i : ±a/3] 
and the two nodes, [0 : ±1 : 1], each with multiplicity 2. We divide these points into 
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two conjugate sets S U S where 

S = {[0 : 1 : 1], [0 : -1 : 1], [2 : : i], [2 : -v^ : i], [2 : i : y/3], [2 : i : -v^]} . 

The vector space of cubics in C [x,y,z] vanishing on these six points is four dimen¬ 
sional and we extend an to a basis {an, (I 12 , «i 3 , Om} for this space, where 

a i2 = ix 3 + 4 ixy 2 — 4 x 2 z — 4y 2 z + 4z 3 , 
a 13 = —3 ix 3 + 4 x 2 y + 4ixy 2 — 4 y 3 + 4 yz 2 , 
an = —x 3 — 2 ix 2 y — 2 ix 2 z + 4 xyz. 

Let M = xA+yB + zC . The two 4x4 polynomial matrix equations aM = (/, 0) and 
Ma T = (/, 0) T give us 120 affine linear equations in the 48 variables Aij, 

For example, the first entry of the vector aM is 

(An + *^2i — “4iA 3 i — An)x^ + ( 4 j 43 i — 2iAn + Bn + il?2i — — Bn)x 3 y 

+ (-2An + 4iAn + 4iA 3 i + 4 B 3 i - 2iBn)x 2 y 2 + (-4T 3 i - 2Bn + 4iB 2 \ + 4iB 3 i)xy 3 
— 41?3iy 4 + (—4^21 — 2 ivl 4 i + CVi + iC 2 i — SiC 3 1 — Cn)x 3 z 

+ (4^41 — 4 H 21 — 2iBn + 4 C 31 — 2iCn)x 2 yz + (—4^21 + 4Bn — 2Cn + 4iC*2i + 4iC 3 \)xy 2 
+ (— 4B 2 \ — 4C 3 i)y 3 z + (—2Tn — 4C*2i — 2iC*4i)x 2 2 2 + (4y43i — 2Bn + 4Cn)xyz 2 
+ ( 4 S 31 — 4C 2 i)y 2 z 2 + ( 44 I 21 — 2Cn)xz 3 + (4i?2i + 4C 3 \)yz 3 + 4 C 212 ; 4 . 


Identifying this polynomial with / gives us 15 affine linear equations. For example, 
from the monomial x 4 , we see that An + ^21 — 3iA 3 i — A 41 = 1. Similarly, from 
each of the other entries of aM — (/, 0) and Ma T — (/, 0) T we get 15 affine linear 
equations in the 3 ■ 4 2 variables A tl , Bij, Cij, for a total of 2 • 4 • 15 = 120. The unique 
solution to these 120 equations gives the Hermitian matrix representation 


xA + yB + zC 


/ 14x 2z 

1 2 z x 

8 —2 ix — 2 y 0 

\— 2 i(y — z) ix + 2 y 


2 ix — 2 y 2 i(y — z)\ 
0 — ix + 2 y 

x ix — 2 z 
—ix — 2 z 4x J 


whose determinant is (1/256) ■ /. 


Ideally we would like to carry out our algorithm symbolically, as in this example, 
but the required field extensions generally will be too large. Given a hyperbolic 
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polynomial / G Q [x,y,z\d, one can ask: What is the field extension necessary to 
carry out the above construction symbolically? 

In fact, after computing the points V(f,g ) and splitting them as S U S, all of 
the remaining steps in the algorithm only require linear algebra, and thus can be 
done with rational arithmetic. Unfortunately, for generic g G M[x, y, z]d- 1 the Galois 
group of the intersection V(/, g) is the full symmetric group of size d(d — 1). It seems 
very hard to say anything about the smallest field extension of the points V(f,g ) 
for special choices of the interlacing polynomial g. 

Remark 5. Since the primary computational task is to find the coordinates of the 
points in V(f,g), it would be worth investigating whether there is a way to find 
the polynomials a±j without computing the points S explicitly. The analogous 1 - 
dimensional problem is, given a positive real polynomial f(t) of degree 2 d , to write 
f(t) as g(t)g(t ) for some degree d polynomial g G C[t]<d. This task is known as 
“polynomial spectral factorization”, see, for example, j2]. Both problems involve 
finding a basis of polynomials (in either C[t]<d or C[x, y, z]d-i/(f)) that vanish on a 
set of non-real points, given a real polynomial that vanishes on both the set and its 
complex conjugate. This higher dimensional analogue does not seem to have been 
considered, but strikes us as worthy of study. 

2. Numerical Implementation 

Here we discuss our numerical implementation of this algorithm in Mathematica 
which is available as supplementary material on the website http://people.math. 
gatech.edu/~rsinn3/numdetrep.html. Below we give computation times and an 
error measure for 100 random tests for polynomials of degrees 3-15. Overall, this 
method results in very fast computations compared to other available methods, 
although the accuracy becomes poor for large d. The numerical accuracy of the ap¬ 
proximation of the common zeros of the hyperbolic polynomial / and the interlacing 
polynomial g (step (Al)) strongly affects the final error in the numerical deter mi - 
nantal representation. The following steps (A4) and (A5) in our algorithm are all 
linear algebra problems that can be solved accurately. 

One issue with numerical computations is that the affine linear equations in (A5) 
are overdetermined - there are d 3 + 3d 2 + 2d equations in 3d 2 variables. With small 
numerical errors, these equations no longer have a solution. In our implementation 
of the algorithm, we solve this by the standard method of taking a least squares 
solution to the system. 

The step that takes the most significant amount of computation time by far is 
computing the points V(f,g). The following steps in our algorithm 

• Computing the basis (an,..., aid)- 

• Translating the polynomial equations (A5) into a system of linear equations. 

• Solving the resulting least squares problem. 

take much less time in our implementation. 

Overall this method finds (approximate) determinantal representations surpris¬ 
ingly fast. To test our code, we generated hyperbolic polynomials of degree d by 
taking the determinant of xl + y(B + B T ) + z(C + C T ) where B and C are random 
d x d matrices whose entries are chosen randomly from a normal distribution with 
mean 1 and standard deviation 0.5. Any such determinant will be hyperbolic with 
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degree 

3 

4 

5 

6 

7 

8 

9 

time (sec) 

0.05 

0.11 

0.25 

0.56 

1.19 

2.36 

4.47 

time V(f,g) 

0.04 

0.09 

0.20 

0.45 

0.96 

1.89 

3.58 

error 

5 • 10- 13 

3- 10" 12 

bO 

h- 1 

O 

1 

t— 1 

bO 

o 

1 

o 

4^ 

1— 1 

o 

1 

cr> 

1 — 1 

o 

1 

00 

o 

1 

si 

relative error 

2 • HT 14 

9 • 10 -14 

2 • 10- 13 

o 

1 

CO 

2 • 10' 12 

h- 1 

H- 1 

O 

1 

9 • 10- 12 


degree 

10 

11 

12 

13 

14 

15 

time (sec) 

9.04 

13.81 

22.95 

37.06 

58.08 

89.71 

time V(f,g) 

6.49 

11.20 

18.72 

30.37 

47.84 

74.31 

error 

2•10“ 6 

0.03 

6 • 10” 4 

0.01 

5.89 

47827 

relative error 

6 • 10” 11 

3•10“ 7 

5 • 10" 10 

4•10~ 9 

5•10" 7 

I— 1 

o 

1 


Table 1 . Average times and errors for computing determinantal representations. 


respect to the point [1:0:0]. Averaging test times for 100 examples in each degree 
gave the computation times shown in Table 1. 

Here “error” means the maximum over the absolute values of the coefficients of 
the difference between the original polynomial / and the appropriately scaled deter¬ 
minant c • det(M). We also found it useful to look at the “relative error”, meaning 
the error divided by the largest coefficient of /. This is a more reasonable measure 
for the accuracy in the coefficients of a hyperbolic polynomial because they tend to 
become very large in higher degree. Indeed, differences between the coefficients of a 
hyperbolic polynomial can grow exponentially with the degree, as in [TJ §4], 

One additional source of numerical errors is the computation of the determinant 
of the output of our algorithm. Because of the size of this matrix, a symbolic compu¬ 
tation of the determinant is infeasible and instead we compute it by interpolation on 
a random set of points on the unit circle. Then we use the interpolated polynomial 
to compute the errors in the coefficients. Numerical testing suggests that numerical 
error in the coefficients of the interpolated determinant are small. 

For comparison, the only other known methods for computing def ini te determi¬ 
nantal representations are discussed in [5]. Here, finding definite determinantal 
representations is already extremely time consuming for quintics (d = 5) and practi¬ 
cally infeasible for larger degrees (d > 6). Thus the method described above provides 
a great improvement in computation ability. 

Remark 6. As we use a random matrix of linear forms to generate our test polyno¬ 
mials, one might also think to compare the starting matrices to those in the output, 
rather then the coefficients of the determinants. The reason this cannot be done is 
that a given hyperbolic polynomial / has an infinite collection of Hermitian deter¬ 
minantal representations / = det(M), even when considered up to the equivalence 
M ~ U*MU for U G GL(d, C). There is no reason for the starting determinantal 
representation xl + y(B + B T ) + z(C + C T ) and the output xMi + yM 2 + zM 3 to 
be equivalent, and so we cannot meaningfully compare their entries. 

Remark 7. The algorithm presented here could also be used to construct a Hermit¬ 
ian determinantal representation of a non-hyperbolic polynomial /, under certain 
conditions. If / G M[x, y, z\d, then one would need a polynomial g G M[x, y, z]d-i 
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such that the intersection of the two curves Vc (/) D Vc (g) is transversal and con¬ 
tains no real points. Furthermore, the vanishing ideal of S in step (A4) of our 
algorithm must have dimension d in degree d— 1; this translates into a condition on 
the corresponding divisor on Vc(/) (see P, §4.5]). It is a non-trivial fact that this 
condition is always met in our original setup when / is hyperbolic and g interlaces 
/. This follows from P Thm. 5.4] or the main result of P; a fuller discussion of 
the non-degeneracy condition in terms of points on the Jacobian can be found in 
p. In the non-hyperbolic case, the condition is at least met for a generic choice of 
g e M[x, y, z]d-i- If d is odd, then the reality condition can be obtained by picking 
g to have no real points, for example by taking g as a sum of squares. Otherwise it 
is not clear how to pick such a polynomial g. 

References 

[1] Alfred Cardew Dixon. Note on the reduction of a ternary quantic to a symmetrical determinant. 
Math. Proc. Cambridge Phil. Society, 11:350-351, 1902. 

[2] Tim N. T. Goodman, Charles A. Micchelli, Giuseppe Rodriguez, and Sebastiano Seatzu. Spec¬ 
tral factorization of Laurent polynomials. Adv. Comput. Math., 7(4):429-454, 1997. 

[3] J. William Helton and Victor Vinnikov. Linear matrix inequality representation of sets. Comm. 
Pure Appl. Math,, 60(5):654-674, 2007. 

[4] Daniel Plaumann, Bernd Sturmfels, and Cynthia Vinzant. Quartic curves and their bitangents. 
J. Symbolic Comput., 46(6):712 733, 2011. 

[5] Daniel Plaumann, Bernd Sturmfels, and Cynthia Vinzant. Computing linear matrix repre¬ 
sentations of Helton-Vinnikov curves. In Mathematical methods in systems, optimization, and 
control, volume 222 of Oper. Theory Adv. Appl., pages 259-277. Birkhauser/Springer Basel AG, 
Basel, 2012. 

[6] Daniel Plaumann and Cynthia Vinzant. Determinantal representations of hyperbolic plane 
curves: an elementary approach. J. Symbolic Comput., 57:48-60, 2013. 

[7] David E. Speyer. Horn’s problem, Vinnikov curves, and the hive cone. Duke Math. J., 
127(3):395-427, 2005. 

[8] Victor Vinnikov. Selfadjoint determinantal representations of real plane curves. Math. Ann., 
296(3) :453-479, 1993. 

[9] Victor Vinnikov. LMI representations of convex semialgebraic sets and determinantal repre¬ 
sentations of algebraic hypersurfaces: past, present, and future. In Mathematical methods in 
systems, optimization, and control, volume 222 of Oper. Theory Adv. Appl., pages 325-349. 
Birkhauser/Springer Basel AG, Basel, 2012. 

Universitat Konstanz, Germany 

E-mail address: Daniel.Plaumann@uni-konstanz.de 

Georgia Institute of Technology, Atlanta, GA, USA 
E-mail address: rsinn3@math.gatech.edu 

University of Michigan, Ann Arbor, MI, USA 
E-mail address: speyer@umich.edu 

North Carolina State University, Raleigh, NC, USA 
E-mail address: clvinzan@ncsu.edu 


